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The objective of this work can be divided into two parts. The first one is to propose an extension of the force density method 
(FDM)|2], a form-finding method for prestressed cable-net structures. The second one is to present a review of various form- 
finding methods for tension structures, in the relation with the extended FDM. 

In the first part, it is pointed out that the original FDM become useless when it is applied to the prestressed structures that consist 
of combinations of both tension and compression members, while the FDM is usually advantageous in form-finding analysis of 
cable-nets. To eliminate the limitation, a functional whose stationary problem simply represents the FDM is firstly proposed. 
Additionally, the existence of a variational principle in the FDM is also indicated. Then, the FDM is extensively redefined by 
generalizing the formulation of the functional. As the result, the generalized functionals enable us to find the forms of tension 
structures that consist of combinations of both tension and compression members, such as tensegrities and suspended membranes 
with compression struts. 

In the second part, it is indicated the important role of three expressions used by the description of the extended FDM, such as 
stationary problems of functionals, the principle of virtual work and stationary conditions using V symbol. They can be commonly 
found in general problems of statics, whereas the original FDM only provides a particular form of equilibrium equation. Then, 
to demonstrate the advantage of such expressions, various form-finding methods are reviewed and compared. As the result, the 
common features and the diff'erences over various form-finding methods are examined. Finally, to give an overview of the reviewed 
methods, the corresponding expressions are shown in the form of three tables. 
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'1. Introduction 

This is a revised version of fH] . 

The objective of the first half of this work is to propose an 
extension of the force density method (FDM)|2], a form-finding 
jnethod for prestressed cable-net structures. Particularly, for the 
prestressed tension structures, form-finding is a process to en- 
sure them to have a prestress state, because the existence of a 
prestress state highly depends on the form of the tension struc- 
ture . 

In section 2, the original FDM is described with its major 
advantage in form-finding process of cable-net structures. In 
addition, it is pointed out that the FDM become useless when it 
is applied to the prestressed structures that consist of combina- 
tions of both tension and compression members, e.g. tensegri- 
ties. Therefore, the FDM has a scope for extension. 

In section 3, a functional whose stationary problem sim- 
ply represents the original FDM is firstly proposed. Addition- 
ally, the existence of a variational principle in the FDM is also 
indicated, although the formulations provided by the original 
FDM look diff'erent from those related to the variational princi- 
ple. The clarified functional enables an extension of the FDM. 
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In section 4, the FDM is extensively redefined by general- 
izing the formulation of the functional. As the result, the gener- 
alized functionals enable us to find the forms of tension struc- 
tures that consist of combinations of both tension and compres- 
sion members, such as tensegrities and suspended membranes 
with compression struts. Moreover, it is pointed out that various 
functionals can be selected for the purpose of form-finding. 

In section 5, some numerical examples of the extended FDM 
are illustrated to show that the newly introduced functionals 
enable us to find the forms of tension structures that consist 
of combinations of both tension and compression members, 
such as tensegrities and suspended membranes with compres- 
sion struts. 

In section 6, in which the second half of this work is de- 
scribed, it is firstly indicated the important role of three ex- 
pressions used by the description of the extended FDM, such 
as stationary problems of functionals, the principle of virtual 
work and stationary conditions using V symbol. They can be 
commonly found in general problems of statics, while the orig- 
inal FDM only provides a particular form of equilibrium equa- 
tion. Then, to demonstrate the advantage of such expressions, 
various form-finding methods are reviewed and compared. As 
the result, the common features and the diff'erences over vari- 
ous form-finding methods can be examined. Finally, to give an 
overview of the reviewed methods, the expressions correspond- 
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ing to them are shown in the form of three tables. 

2. Force Density Method 

2.1. Original Formulation 

The FDM is one of the form-finding methods for cable- 
net structures which was first proposed by H. J. Schek and K. 
Linkwitz in 1973. When it is explained, two unique points are 
usually pointed. The first one is the definition of the force den- 
sity and the second one is the linear form of the equilibrium 
equation provided by the FDM. 

As the first one, the force density qj is defined by 

^j = f^j/Lj, (2.1) 

where nj and Lj denote the tension and length of the j-th mem- 
ber of a structure respectively, as shown in Fig. 12. H a). In the 
FDM, each tension member is assigned a positive force density 
as a prescribed parameter, even though nj and Lj are unknown. 
However, In Ref. f?], there is no mention of method to deter- 
mine them. Then, it is sometimes pointed out that some trials 
must be carried out to obtain an appropriate set of force densi- 
ties. 

As the second one, although the form-finding problems usu- 
ally formulated as a non-linear problem, the self-equilibrium 
equation provided by the FDM is formulated as a set of si- 
multaneous linear equations. In detail, when the force densi- 
ties and the coordinates of the fixed nodes are prescribed, the 
self-equilibrium equation of a cable-net structure is expressed 
as follows: 
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where D is the equilibrium matrix and x, j, and z are the col- 
umn vectors containing the coordinates of the nodes. The terms 
with the subscript / refer to the fixed nodes, whereas those with 
no subscript are for the free nodes. 

Using the inverse matrix of Z>, the nodal coordinates of the 
free nodes can be simply obtained as follows: 
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because, in Eq. (12.21) , only jc, j, and z contain the unknown 
variables. 

Because Eq. (12.3b simply represents the common procedure 
to solve a set of simultaneous linear equations, the FDM can be 
easily implemented by general numerical environments. This 
can be a major advantage in form-finding analysis of cable-net 
structures. 

Once the nodal coordinates are obtained, the tension in each 
cable is calculated by using Eq. (12.11) . The obtained set of ten- 
sion represents a self-equilibrium state of the form, i.e. 

n = {qiLw- ,qmLm}, (2.4) 



where m denotes the number of the members. Generally, such 
a form is called a self-equilibrium form and can be used as a 
prestressed structure. 

Using the FDM, as shown in Fig. I2.ir b), the form of a 
cable-net can be varied by varying the prescribed coordinates 
of the fixed nodes and the force densities of the cables. 

n 

^ /l 

q = n /L 

(a) Definition of force density 




(b) Fonn-finding analysis using FDM 



Figure 2.1: Force Density Method 



2.2. Limitation of FDM 




(a) (b) 



Figure 2.2: X-Tensegrities 

In this subsection, the limitation of the FDM is discussed. 
When it is applied to self-equilibrium systems that consist of 
a combination of both tension and compression members, e.g. 
tensegrities, some difficulties arise. 

In detail, although it seems possible to assign negative force 
densities to the compression members and positive force densi- 
ties to the tension members, the FDM can not keep its concise- 
ness any longer as discussed below. 

Let us consider form-finding of a prestressed structure which 
is called X-Tensegrity. Two diff'erent forms of X-Tensegrity are 
shown by Fig. 12.21 (a) and (b). An X-Tensegrity is a planar pre- 
stressed structure that consists of 4 cables (tension) and 2 struts 
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(compression). As in the case of general tensegrities, the cables 
connect the struts and the struts do not touch each other. 

For such self-equilibrium systems, due to the absence of the 
fixed nodes, Eq. (12.2b reduces to a simpler form: 



ing complementary solution are as follows: 



D x = 0, Z) J = 0, D z = 0. 



(2.5) 



When D is a regular matrix, because it is obvious that D ^ 
- 0, only the trivial solution, i.e. 



X = 0, J = 0, z = 0, 



(2.6) 



is obtained, which implies that every nodes meet at one point, 
namely [0, 0, 0]. 

On the other hand, when Z> is a singular matrix, i.e. detZ> = 
0, Eq. (12.2b generally has complementary solution, which states 
the possible forms of the structure. Such solutions are obtained 
by analyzing the null space of D. Various methods to analyze 
the null space of D. Various methods have been proposed to 
analyze such a space (see Ref.|3-6]). However, even if the 
complementary solutions can be obtained by such methods, the 
major advantage of the EDM, that the equilibrium equation can 
be simply solved by inverse matrix, vanishes. 

Let us see a simple example, the form-finding analysis of 
X-Tensegrity which is shown by Eig. 12.21 When the EDM is 
applied to this type of structure, D is calculated by 
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(2.7) 



(2.8) 



(2.9) 



where C is the branch-node matrix (see Ref.[2] for more detail), 
qi," ' , ^4 are the prescribed force densities of the cables, and 
qs, q^ are of the struts. Then D is represented by: 
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(2.10) 



Based on Eq. (12.10b , the detail of the form-finding analysis of 
X-Tensegrity is as follows: 

• When the assigned force densities, qi, - - • ,qe, are in the 
proportion 1:1:1:1:-1:-1, D becomes a singular matrix 
having 3 dimensional null-space. Then, many solutions 
are obtained. The components of D and the correspond- 
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(2.11) 
(2.12) 



where a, - • • , / are arbitrary real numbers. This implies, 
for example, that both Eig. I2.2r a) and (b) satisfy Eq. 
(12.5b . The first terms of the right hand sides denote the 
position of the center point, namely [a, J, g\ and the other 
terms state some symmetries that all the solutions must 
have. Note that the particular solution is just jc = j = z = 
0. 

• When the assigned force densities, ^i, • • • , ^6, are not in 
the proportion 1:1:1:1:-1:-1, D also becomes a singular 
matrix but having only 1 dimensional null- space. Eor ex- 
ample, if the force densities are in the proportion 2:2:2:2:- 
1:-1, the components of D and the corresponding com- 
plementary solution are as follows: 
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,^6} = {2,2,2,2,-1,-1}) 



(2.13) 
(2.14) 
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where (2, /?, c are arbitrary. This implies that all the nodes 
meet at one point, namely [a, b,c]. 



3. Variational Principle in the FDM 

Let us consider a simple functional 

U(x) = Yj^jL](x), 



(3.1) 



where wj and Lj denote an assigned positive weight coefficient 
and a function to give the length of the 7-th tension member, re- 
spectively. The column vector x represents unknown variables. 
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which are x, y, and z coordinates of the free nodes. It is gener- 
aUzed as an unknown variable container by 



(3.2) 



where n denotes the number of the unknown variables. Note 
that the coordinates related to the fixed nodes are eliminated 
from X beforehand and directly substituted in Lj. 

Actually, the FDM can be simply represented by Eq. (13.11 ): 
the reason is as follows. 

Let V be the gradient operator by 



df 



dxi 



' dXn 



(3.3) 



which points the direction of the greatest rate of increase of /. 
Let Sx be an arbitrary column vector by 



6x = 



(3.4) 



which is called the variation of x. Then, the variation of a func- 
tion f(x) is defined by 



Sm^Vf'Sx. 



(3.5) 



Taking the variation of Eq. (13. IK the stationary condition of 
the functional is calculated as follows: 
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(3.6) 
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In particular case that {xi, • • • , x„} represents the Cartesian 
coordinates of the free nodes, each Lj may defined by the fol- 
lowing form: 

Lj{px, Py. Pz^ qz^ qy. qz) 

= ^l(Px - qxY + {Py - qyf + {Pz - qzf^ 



where p, q denote two ends of j-th member and px 
note 6 coordinates chosen from {xi , • • • , 



(3.10) 

, de- 
In this case, VLj 
represents two normalized vectors attached to both ends of j-th 
member, as shown in Fig. B.lT a). 

On the other hand, suppose the same member resisting two 
nodal forces applied to both ends, as shown in Fig. B.lT b). If 
the magnitude of the tension of the member is denoted by nj, 
then the magnitudes of the two nodal forces are also nj. 

By comparing Fig. B.lT a) and (b), a general form of the 
self-equilibrium equation for prestressed cable-net structures is 
obtained as 

J^«,VL,=0. (3.11) 

j 



To obtain another general form, taking the inner product of 
Eq. (13.111) with dx, the Principle of Virtual Work for such 
structures is obtained as 



(3.12) 



where SLj is the variation of Lj. 
When a set of nj, i.e. 



(3.13) 



where m denotes the number of the members, satisfies Eq. (13.111 ), 
such a set of nj represents a self-equilibrium state of the struc- 
ture. 

Remembering the definition of the force density, namely 



(3.14) 



Eq. (l2.1l) , Eq. (13.111) can be rewritten as 



which is an alternative form of equilibrium equation provided 
by the FDM. 

Comparing Eq. (13.91 ) and Eq. (13.141) , when Eq. (13.91 ) is 

considered as a equilibrium equation, wj is just a half of qj. 
Moreover, when Eq. (13.11 ) is stationary with a form, it is also 
the result obtained by the FDM when the prescribed distribution 
of qj is as same as wj. 

Therefore, Eq. (13.11 ), whose stationary condition is Eq. (13.91 ), 
is one of the functionals that simply represents the FDM. In ad- 
dition, it is assumed that the assigned weight coefl&cients would 
play the same role in form-finding analysis as the force densi- 
ties in the FDM. 

Because the left hand side of Eq. (13.14b simply represents 
the gradient of Eq. (13.11) , the stationary problem of Eq. (13.11) 
can be solve by general direct minimization approach, such as 
the steepest decent method or the dynamic relaxation method 

liMl. 

Although Eq. (12.21) and Eq. (13.141) look very diff'erent, they 
are accurately identical when each function Lj is defined by Eq. 
(13.101) . Then, let us examine Eq. (13.141) for further comprehen- 
sion of the linear form of equilibrium equation provided by the 
FDM. If the non-zero components of VLj is split out as 
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the components of VLy are calculated as 
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Here, it can be noticed that Lj(x) makes VLy non-linear. Then, 
a linear form can be obtained by multiplying VLy with Lj(x), 
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which is the foundation of the Hnear form of the equiUbrium 
equation provided by the FDM. 

Let us consider a case that each variable Xi is also a function 
of another set of variables {ji, • • • , j^}, i.e. 



xi = xi(yi 



In this case, the variation of x is given by 



Sx = D Sy, 
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On the other hand, the relation between two types of gradients, 
namely with respect to x and j, is given by 
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which implies that the expressions such as Eq. (13. 9K Eq. (I3.11K 
Eq. (13.121) and Eq. (13.141) remain valid when {xi, • • • , x„} rep- 
resents other coordinate, such as the polar coordinate. 

In this fashion, Eq. (l3.14b is the general form of the equi- 
librium equation provided by the FDM. On the other hand, the 
equilibrium equation provided by the original FDM is one of 
the particular forms of the general form, which is only valid for 
the Cartesian coordinate. 

Taking the inner product of Eq. (13.141 ) with Sx, the Princi- 
ple of Virtual Work for the FDM is obtained as 
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Similarly, the Principle of Virtual Work is also deduced from 
Eq. ([321) as: 

Sw = Yj ^^jLjSLj = 0. (3.26) 



As the result, as well as in the general problems of statics, the 
variational principle for the FDM is simply represented by 
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where SYl is defined by 



= vn • Sx. 



(3.27) 
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To conclude this section, it is important to note that, in the 
original reference I*?], Eq. (13.11 ) have been mentioned by the 
following theorem: 

"THEOREM 1. Each equilibrium state of an un- 
loaded network structure with force densities qj is 
identical with the net, whose sum of squared way 
lengths weighted by qj is minimal. " 




(b) Equilibrium state 
Figure 3.1: Linear Member 

4. Extended Force Density Method 

4.1. Generalized Formulation of Functional 




Figure 4.1: X-Tensegrities and Prescribed Parameters 

In this subsection, the FDM is extended for form-finding 
of structures that consist of combinations of both tension and 
compression members, e.g. tensegrities. 

Let us reconsider the form-finding of X-Tensegrity again. 
Although it seems possible to assign negative weight coeffi- 
cients to the compression members and positive weight coef- 
ficients to the tension members, the same difficulties which is 
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Figure 4.2: Relation between Original FDM and Extended FDM 



pointed out in subsection 2.2 also arise from the stationary prob- 
lem of n(x) = 2 In detail, when the assigned weight 
coefficients wj are in the proportion for the 4 ten- 
sion members and the 2 compression members respectively, the 
stationary points forms a space. On the other hand, when wj are 
not int the proportion l:l:l:l:-l:-l,the stationary points vanish. 

First, it is obvious that, without no constraint conditions, ev- 
ery lengths of the members become simultaneously or infinite. 
This is due to the absence of information about the scale of the 
structure. Remembering that, in the original FDM, such infor- 
mation is given by the prescribed coordinates of the fixed nodes, 
let the lengths of the compression members be prescribed. Then, 
using the LagrangeJ\fs multiplier method, a modified functional 
is obtained as 



where the first sum is taken for all the tension members and the 
second is for all the compression members. In addition, Xk and 
Lk denote the Lagrange^fs multiplier and the prescribed length 
of the ^-th compression member, respectively. Note that the 
positive weight coefficients wj are assigned to only the tension 
members and the prescribed lengths Lk are assigned to only the 
compression members as shown in Fig l4.1l 

However, Eq. (14.11) does not completely eliminate the above 
mentioned difficulties. For example, if the assigned weight co- 
efficients of the tension members wi, • • • , W4 are in the pro- 
portion 1:1:1:1, and the prescribed lengths of the compression 
members L5, are in the proportion 1:1, both Fig. 14.1 f a) and 



(b) satisfy the stationary condition of Eq. (14.11) . By using the 
Pythagorean theorem, i.e. = -\-b^, it can be easily verified 
that the sum of squared lengths of the tension members takes 
the same value for both Fig. K.lT a) and (b). Then, it is assumed 
that such difficulties depend on the power of L^, i.e. 2. 
Thus, other functionals, such as 

n (X, i) = ^ wjL'j (x) + {Lk (x) - 4), (4.2) 



are introduced, because it is possible to use other powers of Lj 
instead of 2. 

Solving the stationary problem of Eq. (14. 2K Fig. 14.11 (a) 
becomes the unique solution when the weight coefficients of 
the tension members wi , • • • , W4 and the prescribed lengths of 
the compression members 1/5,1/6 are in the proportion 1:1:1:1 
and 1:1 respectively. On the other hand, when they are 1:8:8:1 
and 1:1, Fig. 14.11 (b) becomes the unique solution. Actually, 
Fig- I4.ir a) and (b) are the real numerical results obtained by 
solving such problems. 

By the way, let us discuss the following general formulation 
of functional: 

n (X, A) = Yj {Lj (x)) ^Yj"^^ (^^ ^-^^ " ^4 ^^-^^ 

j k 

The stationary condition of Eq. (14.31) with respect to x is as 
follows: 

dU ^^^j{Lj(x)) ^ 



Because Eq. (14.4b has the same form of Eq. (I3.11K it can be 
considered as a equihbrium equation. Then, when Eq. I4.3l is 
stationary, the following non-trivial set of axial forces must sat- 
isfy the general form of equilibrium equation: 



' dL,„ 



9 ^^m+r 

(4.5) 

which represents a self-equilibrium state of structure, where m 
and r denote the numbers of the tenion and the compression 
members respectively. 

On the other hand, the stationary condition of Eq. (14.31 ) with 
respect to X is given by 
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Therefore, any functional that compatible to Eq. (14.31 ) has 
a possibility to be used for such form-finding problems. From 
now on, let us call nj the element functional. Then the following 
policy is proposed: 

• Perform form-finding analysis by solving a stationary prob- 
lem that is formulated by freely selected element func- 
tionals. 

Taking the inner product of Eq. (14.41 ) with Sx, the Principle of 
Virtual Work is obtained as: 



(4.7) 



d7:j{Lj(x)) 
Sw = SL. — + Zj ^^^^^ = ^' 

r ^ k 



Additionally, replacing the partial derivatives in Eq. 14.71 
by Hj, the following form can be also used as the Principle 
of Virtual Work for general prestressed structures that consist 
of combinations of both tension and compression members: 

6w = Yj ^A^j^^^j + Z ^^^^^ " ^^-^^ 

j k 

Comparing Eq. (14.71) and Eq. (14.81) , if w^L^ is selected as 
the element functional, the following relations are derived: 



■' J 
~dL~ 



^^jLj, = nj/2Lj. 



(4.9) 



Hence, wj can be considered as a half of the force density of the 
7-th member. 

On the other hand, if WjL] is selected, then. 



•' J 



■ 4w/L 



fijlAL]. (4.10) 



Thus, in this fashion, various quantities that are similar to the 
force density can be defined. Then, let us call the new quanti- 
ties, such as Wj = fijlAL^j, the extended force density. 

Apartting from the linear form of the equilibrium equation, 
now, the main characteristics of the original EDM are reconsid- 
ered as follows: 



• The coordinates of the fixed nodes are prescribed as con- 
straint conditions. 

• The force densities qj = nj/Lj are assigned to each ten- 
sion member as known parameters. 

On the other hand, for example, when w^Lj is selected as the el- 
ement functional, the main characteristics of the extended EDM 
are as follows: 

• The coordinates of the fixed nodes and the lengths of the 
compression members are prescribed as constraint con- 
ditions. 

• The extended force densities, e.g. wj = nj/4L^j, are as- 
signed to each tension member as known parameters. 

Therefore, the extended EDM can be considered as similar method 
to the original EDM. 

Considering both approach as solving the stationary prob- 
lems, their main diff'erence is related to the form of the station- 
ary conditions and the selection of the computational methods. 
In the original EDM, they are as follows: 

• The stationary condition of functional is represented by a 
particular form. 

• The stationary condition is simply solved by using an in- 
verse matrix D~^. 

On the other hand, in the extended EDM, they are as follows: 

• The stationary condition of functional is represented by a 
general form. 

• The stationary condition is solved by general direct min- 
imization approaches. 

As an overview of the relation between the original and the ex- 
tended EDM, Eig. 14.21 shows a diagram of both procedures. 

4. 2. Additional Analyses 

In this subsection, some additional numerical analyses are re- 



(-10,-10,0) 



(10,-10,0) 



(-10,10,0) (10,10,0) 

Figure 4.3: Analytical Model 

ported to supplement the concept of the extended EDM. 

Let us consider a net that consists of 220 cables (tension 
members) connecting one another and having 5 fixed nodes as 
shown in Eig. 14.31 The prescribed coordinates of the fixed 
nodes are also shown in the figure. 
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(i)ZLj min 



(ii) ^Lj ^ min (iii) 2 ^] ^ rnin 

Figure 4.4: Optimization Results of Cable-nets 



(iv) 2 ^7 ^ min 




(i) TjLj ^ min (ii) X ^ min (iii) YL^. ^ min (iv) X ^ min 

•' J J J 

Figure 4.5: Optimization Results of Simplex Tensegrities 



Next, let us find the forms taking minimum numbers of 
liy • • • ' namely 



5.7. Structures Consisting of Cables and Struts 



z 



J 



min, {1,2,3,4}), 



(4.11) 



where Lj denotes the length of the j-th cable. The results of 
minimization processes are shown in Fig. 14.41 

On the other hand. Fig. 14.51 shows the other results of the 
same series of minimization processes performed on another 
model, which is based on Simplex Tensegrity. A Simplex Tenseg- 
rity is a prestressed structure that consists of 9 cables(tension) 
and 3 struts(compression). In addition, the minimization pro- 
cesses were only performed on the cables, whereas, the lengths 
of the struts were kept constant at prescribed length, 10.0, dur- 
ing the processes. 

Comparing particularly Fig. I4.4f ii) and Fig. I4.5f ii), w^L^ 
seems not good for form-finding of tensegrities. 



For more detail, when L, 



0, VLj can no be defined be- 



cause VL becomes division by zero (see Eq. (13.161) ). Therefore, 
three of the results, i.e. Fig. l4lT i), Fig. l43l i) and Fig. lOii). 
are only the solutions of minimization problems, whereas the 
others are also the solutions of stationary problems. 

5. Numerical Examples 

In this section, numerical examples of the extended FDM 
are presented. 

In the examples, the stationary problems are represented in 
the following form: 

n (X, A) = n^(x) + [U (x) - Lk). (5.1) 

k 

Then, for simplicity, the problems were solved by general direct 
minimization approaches, in which just IVwix) were minimized 
as objective functions and the lengths of the struts were kept 
constant at the prescribed lengths Lk during each minimization 
process. Hence, only x, or the form, was obtained in each prob- 
lem. 





(a) Connection (b) Result 

Figure 5.1: Simplex Tensegrity 



35 



37 39 1 3 



36 



38 



40 2 



Figure 5.2: Sequential Numbers on Struts 

As mentioned in section 4.2, a form of the Simplex Tenseg- 
rity that consists of 9 cables(tension) and 3 struts(compression) 
can be obtained by solving the following problem: 

n (X, a)=Yj^] (•^) + Z (^^ ^-^^ " ^'^) 

j k 

Stationary. (5.2) 



Here, in the relation with Eq. (15.11 ), the objective function 
Hv, is Y,jL]{x). 

The Principle of Virtual Work corresponding to Eq. (15.21 ) 

is as follows: 



dw ■■ 



2 ^L]dLj + ^kSLk = 0. (5.3) 



Cable#(wi) 


Node# 


Cable# (W2) 


Node# 


1 


1-3 


41 


1-4 


2 


2-4 


42 


2-5 










39 


39-1 


79 


39-2 


40 


40-2 


80 


40-3 



(Ci) 



Cable#(wi) 


Node# 


Cable# (W2) 


Node# 


1 


1-5 


41 


1-6 


2 


2-6 


42 


2-7 










39 


39-3 


79 


39-4 


40 


40-4 


80 


40-5 



(C2) 



Cable#(wi) 


Node# 


Cable# (W2) 


Node# 


1 


1-19 


41 


1-20 


2 


2-20 


42 


2-21 










39 


39-17 


79 


39-18 


40 


40-18 


80 


40-19 



(C9) 



Table 1 : Connections by Cables 



In the analysis, every prescribed lengths of the struts, Ly^, 
were set to 10.0. The connection between the struts and the 
cables in a Simplex Tensegrity is as shown in Fig. 15.11 (a). The 
obtained result is shown by Fig 15. 11 (b). 

Generally, in the direct minimization approaches (see Ref. 
different initial configurations of x may give diff'erent 
results, because the functionals are basically multimodal. 

Then, diff'rent random numbers from -2.5 to 2.5 were roughly 
set to the initial configuration of x in each analysis in order to 
obtain local minimums as many as possible, because it is not 
only the global minimum but any local minimum has an ability 
to be used as a tension structure. 

In this example, particularly, only Fig. 15.11 (b) were con- 
stantly obtained. However, the same strategy was used in the 
following examples and in some of them, many local minimums 
were obtained. 

Let us consider more complex tensegrities such as a system 
that consists of 80 cables (tension) and 20 struts (compression). 
Let us assign sequential node numbers to all the ends of the 
struts, as shown in Fig. 15.21 

Even there are a variety of connections between the struts by 
the cables, 9 of connections were tested. For each connection, 
the node numbers that each cable connects are as shown in Tab. 

□ 

In this example, the following stationary problem was for- 



mulated and a series of form-finding analyses were carried out: 

40 80 

n (X, ^1^' ^-^^ + Z ^^^'j ^-^^ + Z (^^ ^-^^ " ^'^) 

7-1 j-41 k 

Stationary, (5.4) 

in which the cables were divided into two groups and wi de- 
notes the common weight coefficients for the first group, whereas 
W2 is for the second group. In addition, every prescribed length 
of the struts, Ly^, were constantly set to 10.0. 

The Principle of Virtual Work corresponding to Eq. (15.41 ) 
is as follows: 



dw ■■ 



40 



^ AwiL]6Lj + Yj ^^2L]SLj + ^ XkSU - 0. (5.5) 



7=1 



7=41 



When wi : W2 = 1 : 2, Fig. 15.31 shows the most frequently 
obtained results for each connection. Fig. 15.41 (j) to (1) shows 
how the form varied when the proportion between wi and W2 
was varied. Interestingly, between Fig. I5.4l (k) and (1), a transi- 
tion of the form was observed. 

It must be noted that the results shown by Fig. I5.3l are just a 
fraction of various obtained results and a lot of local minimums 
were obtained for each connection, which implies that the func- 
tionals are multimodal. An example of such local minimums 
are given by Fig. 15.41 (F) and (m). Although both results have 
exactly the same connection and the prescribed parameters, ex- 
cept the initial configuration of x, their forms look completely 
diff'erent. This is due to the random numbers which were set to 
X in each initial step. 

5.2. Structures Consisting of Cables, Membranes and Struts 
HL-L) 




Figure 5.5: Analytical Model 

For form-finding of structures that consist of combinations 
of cables (tension), membranes (tension), and struts (compres- 
sion), if the cables are represented by a set of linear elements 
and the membranes, by a set of triangular elements, Eq. (14.31 ) 
can be extended as follows: 

n (x, ^)=Yj i^j M) + X ^^^^^ 

j k 

^^/(L/(x)-L/), (5.6) 
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(a)C^,w^:w-l:2 



(b)q, w-w -1:2 (c)C3, w,:>t; -1:2 (d)^, >v^:w =1:2 {q)C^,w^:w -1:2 






(f)q, H',:w,-1:2 (g)C,, w,:w =1:2 (h)C,, h^,:w3-1:2 (i)C,, w,:w =1:2 



Figure 5.3: Discovered Tensegrities 
Effect of Parameter Change ------- 




Local Minimums 





Q)C^,w-w=\:4 {k)C^,w;.w=\'A (1)^, w,:w= 1:0.5 , , (f) C^,w^:w =1:2 (m) C^,w;.w=\:2 

Figure 5.4: Variety of Forms (C^) 



where the first sum is taken for all the linear elements, the sec- 
ond is for all the triangular elements and the third is for all the 
struts. In addition, Lj and Sk are defined as the functions to 
give the length of the j-th linear element and the area of the 
^-th triangular element respectively. 

The stationary condition of Eq. (15.61) with respect to x is as 
follows: 



'dx ~ ^~ ^ 



dLj 



E 



dSk 



V^, + ^^/VL/ = 0. 



k " I 
Replacing the partial diff'erential factors by 

dnj{Lj(x)) dnkiSkix)) 



(5.7) 



a general form that can be considered as a self-equilibrium equa- 
tion for such systems is obtained as: 



In order to alter the cables in the tensegrities by tension 
membranes, a form-finding analysis based on the above formu- 
lations was carried out with an analytical model shown by Fig. 
15.51 The model is based on the cuboctahedron and consists of 
24 cables, 6 membranes, and 6 struts. In detail, every members 
were translated to purely geometric components such as curves, 
surfaces and lines, then, each curve were discretized by 8 lin- 
ear elements and each surface was discretized by 128 triangular 
elements. 

In the analysis, the following stationary problem was for- 
mulated and solved: 

n (X, A)=Yj ^jL'j ix) ^Yj'^kSl 

j k 

^ Ai [Li ix) - Li^ stationary. (5.10) 



The Principle of Virtual Work corresponding to Eq. (15.101 ) is 

as follows: 



Yj ^j^^j Yj ^^^^ ^ Yj ^^^^^ " ^ ^^-^^ = X + Y ^^kSkSSk + Y ^^^^^ = ^- ^^-^ ^) 



Taking the inner product of Eq. (15.81 ) with Sx, the Principle 
of Virtual Work corresponding to Eq. (15.81 ) is obtained as fol- 
lows: 

sw=Y ""j^^j + Y "^^^^^ + Y ^^^^^ " ^- ^^-^^ 



At first, all of the weight coefficients of the linear elements 
were set to 2.0, those of the triangular elements, 1.0, and the 
prescribed lengths of the struts, 10.0. Then the initial result 
shown by Fig. I5.6l (n) was obtained. By varying wy, Wk and L/, 
the form was able to be varied as shown in Fig. l5.6l (o) to (q). 
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(n) Initial Result 



(o) Revision 1 



(p) Revision 2 



(q) Final Result 



Figure 5.6: Result of Complex Tensegrity 





(r) Initial Result 



(s) Parameter Change 1 (t) Parameter Change 2 (u) Parameter Change 3 (v)Parameter Change 4 (w) Final Result 



Figure 5.7: Form-finding of Suspended Membrane 



5.3. Structures Consisting of Cables, Membranes, Struts and 
Fixed Nodes 

A form-finding analysis of a suspended membrane struc- 
ture based on the famous Tanzbrunnen was carried out. It is 
located in Cologne (Koln), Germany, and was designed by F. 
Otto (1957). 

In the analysis, the following problem was formulated and 
solved: 

n (X, = ^i^l (•^) + Z "^^^^ ^-^^^ 

j k 

^ Ai [Li (x) - L/) stationary, (5.12) 

/ 

where, as well as in the previous example, the first sum is taken 
for all the linear elements, the second is for all the triangular 
elements, and the third is for all the struts. As well as in section 
3, the prescribed coordinates of the fixed nodes are eliminated 
from X beforehand and directly substituted in Lj and S k- 

By varying wj, Wk and L/, as shown in Fig. 15.71 the form 
was able to be varied. Note that Fig. I5.7r w) looks having a 
close form to the real one. 

6. Review of Various Form-Finding Methods 

In the description of the extended FDM, which is just intro- 
duced in the previous sections, three diff'rent types of expres- 
sions are mainly used, they are, stationary problems of func- 
tionals, the principle of virtual work, and stationary conditions 
using V symbol. Such expressions can be commonly found in 
general problems of statics. 

In this section, by using such expressions, various form- 
finding methods are reviewed and compared, in the relation 
with the extended FDM. The methods to be reviewed are, the 
original FDM, the surface stress density method (SSDM) [7], 
and the methods to solve the minimal surface problem, a varia- 
tional method for tensegrities fs*] . 

First, let us review the SSDM, which is also an extension of 
the FDM and a form-finding method for membrane structures. 
It was proposed by B. Maurin et al., in 1998. 



In the SSDM, the membranes are discretized by many tri- 
angular membrane elements and in each elements, the Cauchy 
stress tensor cr^ is assumed as uniform and isotropic, i.e. cr^ = 
(T^^, in order to obtain uniform stress surfaces. As an analogy 
of the definition of the force density, the surface stress density 
Qj in each element j is defined by 



(6.1) 



where ctj is just the scalar multiple of &j with the element thick- 
ness tj and S j denotes each element area. Then, an equilib- 
rium equation is formulated by considering the equilibrium of 
all nodes of the triangular elements. 

Let us rewrite the equilibrium equation provided by the SSDM 
by using V symbol, which is the same fashion that applied to 
the original FDM (see section 3). First, let S(x) be a function 
to give the area of a triangle determined by three nodes whose 9 

coordinates are included in x = |^ xi • • • x„ j . When VS 
is defined by 

dS dS 

dxi ' ' dxn 



VS = 



(6.2) 



it represents three vectors attached to each node, as shown in 
Fig.EIta). 

By the way, let a triangular membrane element, of which the 
thickness is assumed as uniform and denoted by t, be resisting 
three nodal forces applied to each node. For the Cauchy stress 
filed in each element, in the same fashion of the SSDM, let 
(T^ = &S^^ and (T = &t. When such an element is in equilibrium 
with the three nodal forces, the nodal forces can be calculated 
uniquely, and it is as shown in Fig. 16.11 (b). 

Comparing Fig. 16.11 (a) and (b), a general form of self- 
equilibrium equation for general systems that consist of such 
elements is obtained as 

J]cr,V5,. = 0, (6.3) 

j 

taking the inner product of Eq. (16.31 ) with Sx, the Principle of 
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Virtual Work for such a system is obtained as: 

Sw = Yj(^j^^j = ^- (6.4) 

j 

By the way, in the SSDM, the surface stress density Qj is 
defined by 

Qj = cTj/Sj, (6.5) 

then substituting Eq. (16.51) to Eq. (16.31) , a general form for the 
self equilibrium equation of the SSDM is obtained as: 

2 QjSjVSj = 0. (6.6) 

j 

Then, one of the functionals that simply represents the SSDM 
is as follows: 

U(x) = Y,^jSj(x), (6.7) 

j 

because the stationary condition of Eq. (16.71) is given by 

^=Vn = Y,2wjSjVSj = 0, (6.8) 

j 

and when Eq. (16.81) is considered as one of the equilibrium 
equations given by Eq. (16. 3K wj can be considered as just a 
half of Qj. In addition, each wj also represents an extended 
force density such as wj = crj/2S j. 

Based on the proposed functionals for the original EDM 
and the SSDM, namely 2; ^^jLj and "^jSj, the SSDM looks 
truly an extension of the original EDM. 

Moreover, based on the corresponding Principle of Virtual 
Works, i.e. 

Sw = Yj ^^jLjSLj = 0, (6.9) 

j 

j 

2wjLj and 2wjS j can be considered as general forces which act 
within the members or the elements and tend to produce small 
change of Lj and S j, respectively. 

In addition, if the Principle of Virtual Works are written 
in the following forms: 

j 

6w = YjWjS[sj) = 0, (6.12) 

j 

then, the extended force densities, wj = njllLj and = crj/2S j, 
can be considered as general forces which act within the mem- 
bers or the elements and tend to produce small change of Lj and 

Next, let us compare the following two problems: 



R R 

— — a 




(a) VS (b) Equilibrium State 



Figure 6.1: Triangular Element 



n(jc) = ^ Sj(x) stationary, (6.13) 

j 

U(x) = Y^Sj(x) stationary, (6.14) 

j 

because, for the minimal surface problem, Y^j^ j is often used, 
whereas, TjjSj simply represents the SSDM when the distribu- 
tion of the surface stress densities is given as uniform. 

By applying both problems to the same numerical model 
shown in Eig. 16.21 2 pairs of results were obtained as shown 
in Eig. 16.31 In addition, such forms are easily observed by a 
soap-film experiment. 

Eirst of all, due to the fact that they are diff'erent function- 
als, it is not obvious that the stationary points given by Eq. 
(16.14b are minimal surfaces. However, the forms of (a-1) and 
(b-1) look identical with (a-2) and (b-2). On the other hand, 
their mesh distributions look dissimilar, i.e. the results given 
hy YjjSj seem to have better mesh distributions in comparison 
with those by 2; ^ j- 

Then, let us see the Principle of Virtual Works, i.e. 

Sw = Yj^^j = ^^ (6-15) 

j 

^w = ^2^/^y = 0. (6.16) 

j 

Then, it can be noticed that, in Eq. (16.16b , the general forces 
which tend to produce small change of S j are proportional to 
S j, which implies that each element is hard to have bigger or 
smaller area compared to the surrounding elements (see Eig. 
16.41) . On the other hand, in Eq. (I6.15K whatever element area 
that each element has, the coefficients of SS j remain always 1. 
Therefore, as long as the total element area is minimum, each 
element is able to have bigger or smaller area compared to the 
surrounding elements. Thus, the diff'erence appeared in Eig. 16.31 
can be well explained by the principle of virtual works. 

The SSDM has been proposed for structures that consist of 
combinations of membranes and cables. When the SSDM is 
applied to such structures, as same as in section 5.2, the ca- 
bles are represented by linear elements and the membranes are 
represented by triangular elements. Then, the force densities 
are assigned to the linear elements and the surface stress densi- 
ties are assigned to the triangular elements. In such cases, the 
SSDM can be simply represented by 
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Figure 6.2: Form-finding Problem of Simple Membrane 




(a-l)2 S h:4.0 area:122 (a-2)2 h:4.0 area:122 




(b-l)i:^ h:6.5 area:186 (b-2)E^^ h:6.5 area:186 

Figure 6.3: Comparison of^S and 




Figure 6.4: Stress State in SSDM 



U(x) = ^ WjLjix) + ^ W;t5'^(x) stationary, (6.17) 

j k 

where the first sum is taken for all the linear elements and the 
second is for all the triangular elements. Fig. 16.51 shows one 
of the results given by solving Eq. (16.171) . The corresponding 
Principle of Virtual Work is as follows: 

Sw = ^ 2wjLjSLj + ^ 2wkS kSS k = 0, (6.18) 
and the stationary condition is obtained as: 

C\T~f 

— = ^ 2wjLjVLj + 2wkSkVSk = 0. (6.19) 

j k 



Eq. (I6.17I )- (I6.19I ) are just simple compositions of corresponding 
expressions related to the original EDM and the SSDM, which 
imply the potential ability of such expressions for extension. 




Figure 6.5: Membrane with Cables 

Next, let us review form-finding methods which have been 
proposed to determin the forms of tensegrities. Particularly, let 
us examine the following two problems: 

U(x,A)=J]^kj{Lj(x)-Ljf 

j 

+ ^ Ak(Lk(x) - Lk) stationary, (6.20) 

k 

U(x,A)=Yj''jL'j(x) 

j 

+ ^ Ak(Lk(x) - Lk) stationary, (6.21) 

k 

where the first sum is taken for all the cables and the second is 
for all the struts. 

In Ref. fs^, Eq. (l6.2Qb is proposed for the form-finding of 
tensegrities. In Eq. (I6.2QK kj and Lj represent virtual stiff'ness 
and virtual initial length of the 7-th cable respectively, which 
do not represent real material but define special (soft) material 
for form-finding analysis. Therefore, as discussed below, an 
appropriate set of Lj is needed. On the other hand, represents 
just the objective length of the ^-th strut. Eig. I6.6r a) shows 
an example of tensegrities which was obtained by solving Eq. 
(l6.2Qb by the authors. 

On the other hand, Eq. (16.211) is one of the stationary prob- 
lems which was just proposed in this work. Eig. I6.6r b) shows 
an example of tensegrities given by solving Eq. (16.211) . 

With respect to the the second sums for the struts, there look 
no diff'erence. 

On the other hand, with respect to the first sums, which are 
for the cables, some diff'erences can be recognized. They are, 
the powers and the terms that are powered. In addition, while 
the first sum of Eq. (16.201) looks an analogy of elastic energy of 
Hook's spring, the first sum of Eq. (16.211 ) looks diff'erent. 

Then, let us see the Principle of Virtual Works correspond- 
ing to Eq. (|6^ and Eq. (|62T]) , i.e. 

dw^Y^ kj {Lj - Lj) SLj -h ^kSLk = 0, (6.22) 

j k 

6w = Yj ^^jL]SLj -h ^kSLk = 0. (6.23) 

k 
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Thus, it can be noticed that, in Eq. ( 16. 221 ). the general forces 7. Conclusions 

kj (Lj - Li) which tend to produce small change of L, are pro- 

^ ' _ / _ ^ In the first part of this work, the extended force density 

portional to (Lj - Lj). Due to the fact that [Lj - Lj) can take method was proposed. It enables us to carry out form-finding 

negative numbers, some of the cables may become compres- prestressed structures that consist of combinations of both 

sion. Then, it can be noticed that an appropriate set of Lj is jg^^j^^ compression members. 

needed to ensure every {Lj - Lj) be positive. The existence of a variational principle in the FDM was 
_ For thispurpose, one of the simplest ideas to determine each ^^^^^^^ ^ functional that simply represents the FDM 
Lj m Eq. dSp for the cables is to setevery Lj as 0_However, proposed. Then, the FDM was extensively redefined by 
when every Lj are set to m Eq. or Eq. (1623, some generalizing the formulation of the functional. Additionally, it 
difficulties arise as mentioned in section 4. Then, to eliminate indicated that various functional can be selected for form- 
the difficulties, one_ of the simplest ideas is to alter the power g^^j^^g j^^^j^^ structures. Then, some form finding analyses 
of the term [Lj - lj) to other numbers such as 4. Thus, the of diff-grent types of tension structures were illustrated to show 
equations used in the extended FDM, such as Eq. and the potential ability of the extended FDM. 
Eq. ( |6.23 |i, emerge. jj^ tjjg second part, various form-finding methods were re- 
in addition, the Principle of Virtual Work corresponding ^-^^^^^ compared in the relation with the extended FDM. 
to Eq. (EIB is also represented in the following form: ^^^^^ three types of expressions such as the principle of vir- 

^^ tual work, which can be commonly found in general problems 

"^i^ + Zj '^^^^^ " °' of statics and are also used in the description of the extended 

^ * FDM, the common features and differences over different form- 

which states that the extended force densities, i.e. Wj = UjlAL], ^"ding methods were examined, 
can be considered as general forces which act within the cables 

and tend to produce small change of Lj. Acknowledgments 

As a result of above discussion, a common feature which 
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j 

Sw = Yj^^j = ^^ (6.27) 

j 

it can be noticed that the general forces which act within the 
elements or the members remain always positive. 




(a) Result by Eq. fOol (b) Result by Eq. ( fOTl 



Figure 6.6: Form-Finding of Tensegrities 

Finally, the stationary problems of functionals, the principle 
of virtual works and the stationary conditions using V symbol, 
which were just compared in this section, are shown in Tab. |2] 
to |4] as an overview. By using those three expressions that are 
usually found in various problems of statics, the common fea- 
tures and the diff'erences over various form-finding methods can 
be examined, as discussed in this section. Moreover, they also 
enable us to combine or extend the methods in natural ways. 
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Table 2: Table of Stationary Problems 



Force Density Method [ 2] 


Minimal Surface e.g. 1 9] 


cables 


membrane 


(a) n = 2 ~^ stationary. 


(b) U = YjS ^ stationary. 


Surface stress density methodlJ] 


membrane 


membrane with cables 


(c) U = YjWS^^ stationary. 


(d) U = J]wL^ + J]wS^^ stationary. 


Variational Method for Tensegrities| 8] 


Extended Force Density Method (Proposed) 


cables and struts 


cables and struts 


(e) n = 2 5 k(L - 1)2 + X A(L - L) ^ stationary. 


(f) n = 2 wL^ + X A(L - L) ^ stationary. 



Table 3: Table of Principle of Virtual Works 



Force Density Method| 2] 


Minimal Surface e.g. 1 9] 


cables 


membrane 


(a) Sw = Y. IwLdL = 


(b) Sw = Y.SS =0 


Surface stress density method| 7] 


membrane 


membrane with cables 


(c) 6w^J]2wS6S^0 


(d) 6w^J] 2wLSL + Y.2wS6S ^0 


Variational Method for Tensegrities| 8] 


Extended Force Density Method (Proposed) 


cables and struts 


cables and struts 


(e) 5w = 2 kLSL + J]A6L^0 


(f) 6w^J] 4wL^6L + J]A6L^0 


Table 4: Table of Stationary Conditions 


Force Density Method| 2] 


Minimal Surface e.g. 1 9] 


cables 


membrane 


(a) f =i:2wLVL = 


(b) f = X V5 = 


Surface stress density method| 7] 


membrane 


membrane with cables 


(c) ^^X2wSWS^0 


(d) ^ = 2 2wLVL + 2 2w5 V5 = 


Variational Method for Tensegrities| 8] 


Extended Force Density Method (Proposed) 


cables and struts 


cables and struts 


(e) f ^Z^^VL + S^VL^O 


(f) ^ = 2 4wL3 VL + 2 ^VL = 
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Appendix A. Gradients 

Appendix A.l. Gradient of Linear Element Length 
Let p and q denote two nodes. Let 



The length of the Hne determined by p and q is given by 





Px 




qx 


p = 


Py 








_ Pz _ 




_ qz _ 



(A.l) 



i^[Px, Py, Pz, qx, qy, qz) 
= ^j(Px - qxf- + {Py - qy)^ + {Pz - qzY' 
If the gradient of L is defined by 



[ dL 


dL 


dL 


dL dL dL ] 






dp,' 


dqx ' dqy ' dq, \ 



VL 

its components are as follows: 

^l^lElZll Py-^y Pz-qz qx-Px qy-Py qz-Pz A 

[l'l'l'l'l'lI 
Let us investigate SL, i.e. 

SL^VL- ^/ 
dq 



(A.2) 
(A.3) 

(A.4) 
(A.5) 

(A.6) 



represent the Cartesian coordinates of p and q. 



As shown in Fig. lA.U dp and dq are firstly projected to the line 
determined by p and q, then, dL is measured on the line. 
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Figure A. 2: Variation of Element Area 



Figure A. 1 : Variation of Element Length 



Appendix B. Some Remarks of Surface Area 



Appendix A. 2. Gradient of Triangular Element Area 
Let p, q, and r be three vertices. Let 





Px 




qx 




rx 


p = 


Py 


, q = 


qy 








_ Pz _ 




_ qz . 




_ rz _ 



(A.7) 



denote the Cartesian coordinates of p, q, and r. 

The area of the triangle determined by p, q, and r is given 

by 



(N ^(q-p)x(r-p)). 
If the gradient of S is defined by 



yic: - \ ds_ ds_ ds_ 

= [ dp,' dpy' dp,' 

its components are as follows: 



dS_ ] 

' dr. J ' 



(A.8) 
(A.9) 



(A. 10) 



(r-q)x 



,(p-r)x 



,(q-p)x 



1 


















1 




"II 


1 
















1 



















1 







1 












1 




1 



















1 







II 












1 





(A. 11) 



where n is defined by 



n = 



N_ 
W\' 



Let us investigate ^5" , i.e. 



6S - -n-((r-q)x6p + (p-r)x6q + (q-p)x 6r) . 



(A.12) 



(A. 13) 



With respect to dp, for example, when dp is orthogonal to the 
element, (r - ^)x^/7 becomes orthogonal to fi, then dS vanishes 
(see Fig. IA.2b . On the other hand, when 6pi^ parallel to the op- 
posite side, (r - q)x dp vanishes, then 6S vanishes. Therefore, 
only the component of dp which is parallel to the perpendicular 
line from p to the opposite side can produce dS . In other words, 
6S is measured on the plane determined by p, q, and r. 



Appendix B.l. Minimal Surfaces and Uniform Stress Surfaces 
The surface area of a surface is given by 



a = I da. 

Ja 

Here, da is called area element and defined by 
= ^detgijdeh 



da: 



'dO' 



(B.l) 



(B.2) 



where gij and O^j are the Riemannian Metric and the local 
coordinate on the surface respectively. 

Using Eq. (IB. 21) , the variation of the surface area, Sa, can 
be calculated and the result is as follows: 



Sa=^ Jj^'^gij ^dctgijdO'de\ 



(B.3) 
(B.4) 



where g^-^ is the inverse matrix of gij. 

By the way, on a membrane, the 2nd Piola-Kirchhoff stress 
tensor and the Green-Lagrange strain tensor are defined by 



S = 



- ^^^n/'g. ®gj.E^ \{g,j - -g,j)-g' -gK (B.5) 



V5et|7 



where T]^ are the components of the Cauchy stress tensor. In 
addition, g^, g\ gij are the dual bases and the Riemannian met- 
ric defined on a reference configuration. 

Then, the Principle of Virtual Work for membranes is ex- 
pressed as: 



6w 



= [tS: 

Ja 



SEdsi 



(B.6) 



where da, a are related to the reference configuration, and t de- 
notes the thickness. Eq. (IB. 6b reduces to the following form: 

Sw= f tr.g'^Sgijdsi, (B.7) 

Ja 

which does not depend on the reference configuration. 

Because Eq. (BA\f can be transformed into the following 
form: 



Sa 



Ja 



(B.8) 
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when t and are uniform on the surface and when = &6[j^, 
where & is also uniform, then 



Sw = t&da :. 5w = <^ da = 0, 



(B.9) 



which is a simple demonstration of the essential identity of uni- 
form stress surfaces and minimal surfaces. 

Appendix B.2. Galerkin Method for Minimal Surface 

When the form of a surface is represented by n-independent 
parameters such as x = |^ xi • • • x„ j , an approximation of 

da= \ g^%,,da = (B.IO) 

Ja 

can be obtained by the Galerkin method and it is as follows: 

da = {^Jj'jVgijds^j -Sx = 0, (B.ll) 
where V is the gradient operator defined by 

V/-ff ••• f 1 (B.12) 

and 6x is the variation of x, or, just an arbitrary column vector. 

When the form is discretized by m elements, the integral 
can be divided into m independent integrals. Hence 



6x, 



where j is the index of each element. 

In each element, remembering the relation of 

Ja j Ja \j 

the following transformation is also correct: 
r ^"^Vg^^da = V fdal , 

Ja j Ja \j 



(B.13) 



(B.14) 



(B.15) 



due to the fact that S symbol is originally defined by ^ when 
€ is the assigned one-parameter to represent the change of the 
form. 

Therefore, when S j is defined as a function to give j-th ele- 
ment area, i.e. 

Sj = fdal , (B.16) 

Ja \j 



then 



Sa = <^ 



Yj^^j -Sx = ^Yj^^J = (^•^'7) 



V J 



which is the stationary condition of 

U(x) = YjSj(x). 



(B.18) 
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